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Abstract 

A  new  flow  field  was  designed  to  search  flow  fields  fitting  polymer  electrolyte  membrane  fuel  cells  (PEMFCs)  better  due  Its  extensible.  There 
are  many  independent  inlets  and  outlets  in  the  new  flow  field.  The  new  flow  field  we  named  NINO  can  extend  to  be  more  general  when  pressures 
at  the  inlet  and  outlet  vary  and  some  usual  flow  fields  will  be  obtained.  A  new  mathematical  model  whose  view  angle  is  obverse  is  used  to  describe 
the  flow  field. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

With  continuously  increasing  need  of  energy  and  concern¬ 
ing  of  environment,  fuel  cell  as  an  effective  and  clean  device  to 
convert  energy  is  attracting  more  and  more  researchers.  Poly¬ 
mer  electrolyte  membrane  fuel  cells  (PEMFCs)  are  considered 
as  highly  promising  power  sources,  especially  for  vehicular 
applications  because  of  its  low-temperature  operation  and  ease 
of  construction.  Precise  design  of  flow  field  and  well-chosen 
operation  condition  can  greatly  improve  the  usage  of  catalyst, 
meliorate  the  wetting  state  and  simplify  PEMFC  systems  [1], 
Parallel  and  interdigitated  flow  fields  [2 — 4]  shown  in  Fig.  la  and 
b  have  their  own  characteristics.  Parallel  flow  field  exhibits  lower 
pressure  difference  and  inhomogeneous  reactant  gas  distribu¬ 
tion  can  easily  occur.  Products  of  the  electrochemical  reactions 
like  water  and  carbon  dioxide  can  clog  single  channels,  and 
liquid  water  can  flood  the  catalyst  layer.  So  liquid  water  flood¬ 
ing  become  serious  when  the  flux  of  gas  is  small  or  electrical 
current  is  high  in  parallel  flow  fields.  The  forced  flow-through- 
the-electrode  characteristic  makes  the  interdigitated  flow  field 
offset  the  disadvantage  of  the  parallel  due  to  that  the  interdig¬ 
itated  flow  field  could  provide  high  transport  rates  of  reactant 
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and  products  to  and  from  the  inner  catalyst  layer.  The  thickness 
of  diffusion  layer  in  interdigitated  flow  field  is  much  thinner  [3] 
than  that  in  parallel  flow  field  resulting  in  that  the  limiting  cur¬ 
rent  and  maximum  power  of  PEMFC  in  interdigitated  flow  field 
is  higher  than  that  in  parallel  state.  However,  the  high  energy 
lost  in  interdigitated  flow  field  is  due  to  the  high-pressure  drop 
throughout  the  PEMFC,  and  the  velocity  at  the  inlet  and  outlet 
is  extremely  high  compare  to  the  mid  of  the  flow  field  because 
the  fluid  all  crossing  the  flow  field  disperses  at  the  inlet  and 
congregates  at  the  outlet. 

In  order  to  offset  the  shortcoming  of  the  two  flow  fields,  a  new 
flow  field  we  named  NINO  flow  field  which  is  the  abbreviation 
of  n-inlet  and  /(-outlet  flow  field,  shown  in  Fig.  lc,  was  formed 
from  combination  of  the  two.  There  are  many  inlets  and  outlets  in 
the  flow  field,  the  quantity  of  inlet  may  be  not  the  same  as  that  of 
outlet  in  the  new  flow  field.  Every  inlet  and  outlet  is  independent, 
and  the  pressures  at  these  places  are  varied  to  extend  the  new  flow 
field  to  be  more  general.  The  simple  situation  with  symmetric 
boundary  condition  of  the  new  flow  field  was  discussed  here. 

A  new  mathematical  model  is  presented  to  describe  the 
new  flow  field.  The  models  used  by  other  researchers  [1-8] 
were  on  the  whole  or  part  way  of  channel  -a-  gas  diffu¬ 
sion  layer  — +  cathode  -a-  membrane  —>  anode  -a-  gas  diffusion 
layer ->  channel.  However,  the  view  here  was  changed  to  the 
obverse  of  PEMFC,  and  a  two-dimension  model  was  build.  This 
model  can  be  used  to  optimize  the  flow  field  presented  above 
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Nomenclature 

List  of  symbols 

A  area  of  the  flow  field  (m2) 

b  constant  in  Eq.  (16) 

cCh  oxygen  concentration  in  channel  (mol  m-  3) 
cgdl  oxygen  concentration  in  GDL  (mol  m-3) 
c0  mol  concentration  of  oxygen  gas  at  the  interface 
of  GDL  and  catalyst  layer  (mol  m-3) 
c,M-  reference  oxygen  concentration  (mol  m-3) 

Coi  mol  concentration  of  oxygen  gas  at  the  interface 
of  GDL  and  channel  (mol  m-3) 
c02  mol  concentration  of  oxygen  gas  at  the  interface 
of  channel  and  bipolar  plate  (mol  m-3) 
coa  mol  concentration  of  oxygen  gas  in  .—direction 
(mol  m-3) 

cWa  mol  concentration  of  water  in  z-direction 

(mol  m-3) 

C2,  C2',  C2"  symmetry  axes 
c;;1 1  oxygen  concentration  at  inletl  (mol  m-3) 

c"’2  oxygen  concentration  at  inlet2  (mol  m-3) 

cJJ?1  water  concentration  at  inletl  (mol  m-3) 

c“2  water  concentration  at  inlet2  (mol  m-3) 

DCh  diffusion  coefficient  in  channel  (m2  s-1) 

Dgdl  diffusion  coefficient  in  GDL  (m2  s- 1 ) 

Dj  diffusion  coefficient  in  a  nonporous  system 
(m2  s- 1 ) 

Dj,r  reference  oxygen  diffusion  coefficient  (m2  s’  1 ) 
D?  effective  molecular  diffusion  coefficient  of 

species  i  in  the  gas  stream  (m2  s-1) 

De0  effective  molecular  diffusion  coefficient  of  oxy¬ 
gen  in  the  gas  stream  (m2  s-1) 

Of  effective  molecular  diffusion  coefficient  of  water 
in  the  gas  stream  (m2  s’  1 ) 
d  1  high  of  the  channel  (m) 
d2  thickness  of  GDL  (m) 

E  potential  of  fuel  cell  (V ) 

Eq  open  circuit  potential  (V) 

F  Laraday  constant  (96,487  C  mol- 1 ) 
i  species  oxygen  or  water/current  density  (A  m-2) 

i  average  current  of  the  whole  fuel  cell  (A  m-2) 

70  exchange  current  density  of  oxygen  (A  m-2) 

inletl  virtual  inlet  of  first  kind 

inlet2  virtual  inlet  of  second  kind 

;(x,  y)  current  function(A  m-2) 

k  transfer  coefficient 

kp  permeability  of  gas  diffuser  (m2) 

A'pi  permeability  of  gas  diffuser  in  channel  (m2) 

&P2  permeability  of  gas  diffuser  under  shoulder  (m2) 

LI  width  of  ports  (m) 

L2  width  of  shoulders  (m) 

L3  length  of  the  whole  flow  field  (m) 

L4  width  of  the  whole  flow  field  (m) 

m  pressure  changes  of  one  port 

M  molecular  weight  (kg  mol- 1 ) 


Ma  molecular  weight  of  oxygen  (kg  mol- 1 ) 

Mw  molecular  weight  of  water  (kg  mol- 1 ) 

Mn  molecular  weight  of  nitrogen  (kg  mol-1) 
n  port  number  of  one  side 

N  the  total  the  independent  condition  number 
Nc h  oxygen  flux  in  channel  (mol  m-2) 

Agdl  oxygen  flux  in  GDL  (mol  m-2) 
outlet  1  virtual  outlet  of  first  kind 
outlet2  virtual  outlet  of  second  kind 
P  total  gas  pressure  (Pa) 

PQ  atmospheric  pressure  (Pa) 
r  internal  resistance  (Q) 

R  gas  constant  (J  mol- 1  K- 1 ) 
s  symmetry  center 

Schi  source  term  in  equation  of  continuity  in  channel 
(mol  m-2  s-1) 

.S-s|,  i  source  term  in  equation  of  continuity  in  GDL 
(mol  m-2  s-1) 

T  operating  temperature  (K) 

Tt  reference  temperature  (K) 

Mini  velocity  at  inletl  (ms-1) 

Min2  velocity  at  inlet2  (ms-1) 

Mout2  velocity  at  outlet2  (ms-1) 

v  velocity  (ms-1) 

x  A'-direction 

y  v-direction 

Greek  letters 
P  water  content 

e  gas  diffuser  porosity  of  the  GDL 

rj  overpotential  (V) 

T]c  cathode  overpotential  (V) 

X  electro-osmotic  drag  coefficient  of  water 

ji  air  viscosity  (N  m-2  s) 

£  1  dimensionless  parameter  in  Eq.  (41) 

£2  dimensionless  parameter  in  Eq.  (41) 

p  density  of  the  fluid  (kg  m-3) 

a  internal  conductivity  of  the  fuel  cell  (S  m2) 

cr  1 ,  u2,  cr3  symmetry  planes 

Superscripts  and  subscripts 
H+  proton 

i  species 

e  effective 

n  nitrogen 

o  oxygen 

w  water 

r  reference 

ini  inletl 

in2  inlet2 

outl  outlet  1 

out2  outlet2 
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Fig.  1.  Air  flow  fields  for  PEMFC  (a)  parallel,  (b)  interdigitated  and  (c)  NINO. 


and  help  us  to  understand  the  whole  situation  of  all  kinds  of  flow 
fields.  It  is  helpful  to  the  job  of  thinning  [9]  PEMFC  because 
the  difference  of  species  is  small  in  ^-direction  in  thin  flow  field. 
The  cathode  side  of  the  new  flow  field  was  studied,  and  some 
effect  of  base  parameter  was  considered. 

2.  New  flow  field  and  mathematical  model 

2.1.  New  flow  field 

The  new  flow  field  shown  in  Fig.  2  has  five  inlets  and  outlets, 
respectively.  The  quantity  of  the  inlet  and  outlet  can  vary  as  you 
need.  Each  inlet  and  outlet  is  independent  and  can  connect  to  the 
pipeline.  The  inlets  and  outlets  can  exchange  with  or  transform 
to  each  other,  then  the  general  flow  field  will  be  got. 

If  the  flow  field  has  2 n  ports  and  each  port  has  m  pressure 
changes,  the  glancing  number  of  the  pressure  change  can  be 
given  as  follows, 

N  =  m2n  (1) 


where  N  is  the  total  independent  condition  number,  if  the  pres¬ 
sures  at  every  port  are  the  same,  the  fluid  neither  can  come  into 
the  flow  field  nor  can  go  out.  In  this  situation,  the  flow  field 
cannot  work  in  the  course  of  nature,  and  the  number  must  be 
subtracted. 

N  =  m2n  -  m  (2) 

Considering  the  symmetry  of  the  flow  field,  the  number  will 
be  smaller  than  the  one  given  above.  In  the  flow  field,  there  are 
three  symmetry  planes  (al,  a2,  a3),  three  symmetry  axes  (C2, 
C2',  C2")  and  one  symmetry  center(s)  shown  in  Fig.  2.  If  every 
port  is  in  the  same  position  after  operations  of  two  symmetry 
elements,  the  two  symmetry  elements  have  the  same  effect  on 
the  flow  field.  For  example,  the  anticlockwise  order  of  the  port 
is  1,  2,  3,  4,  5,  10,  9,  8,  7,  6  before  operation,  and  became  5, 
4,  3,  2,  1,  6,  7,  8,  9,  10  after  operation  of  al  or  C2.  So  the  al 
and  C2  have  the  same  effects  on  the  NINO,  and  we  only  need  to 
consider  one  of  them.  The  operation  of  a2  and  C2',  C2"  and  s 
have  the  same  condition.  After  the  operation  of  a3,  every  port  is 
in  the  former  position  and  the  result  will  not  be  effected  by  this 
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operation.  As  a  result,  we  only  need  to  consider  the  operation  of 
crl,  <j2  and  s. 

When  the  ports  are  set  to  some  pressures,  after  the  operation 
of  C2,  the  condition  become  superposable  with  that  before  the 
operation,  then  this  condition  only  appears  one  time  under  this 
operation.  If  the  condition  is  not  superposable  with  the  former, 
this  condition  appears  twice  under  this  operation.  The  things 
above  will  happen  in  both  cr2,  C2'  and  C2",  s.  These  conditions 
appearing  repeatedly  must  be  subtracted  from  the  quantity  cal¬ 
culated  from  Eq.  (2).  But  if  the  results  are  the  same  after  the 
operation  of  C2  and  C2',  or  C2  and  C2",  or  C2'  and  C2",  it  is 
needed  to  add  the  number  calculated  here.  Then  the  independent 
condition  can  be  described  as  follows, 

N  =  m2n  -  m  -  N1  -  N2  -  N3  +  N4  +  N5  +  N6  (3) 

where  N 1  is  the  unsuperposition  number  after  operation  of  C2 
or  cj  1 ,  N2  is  the  unsuperposition  number  after  operation  of  C’2' 
or  ct2,  N3  is  the  unsuperposition  number  after  operation  of  C2" 
or  cr3,  N4  is  the  superposition  number  after  operation  of  C2  and 
C2',  N5  is  the  superposition  number  after  operation  of  C2  and 
C2",  N6  is  the  superposition  number  after  operation  of  C2'  and 
C2". 

If  the  pressures  at  the  ports  are  changed,  some  interesting 
results  will  be  got.  First,  if  the  port  6  and  port  7,  2  and  3,  8  and 
9,  4  and  5  are  combined,  respectively,  the  NINO  flow  field  will 
become  the  serpentine  flow  field  [5].  Second  if  the  ports  1,  2,  3, 
4,  5  are  combined  directly,  and  ports  6,  7,  8,  9,  10  are  combined 
too,  but  the  pressure  of  the  two  kinds  of  ports  is  different,  the 
NINO  flow  field  will  be  the  parallel  one.  Third  if  the  ports  1,  3, 
5  have  the  same  pressure,  and  7,  9  have  the  same  pressure,  but 
the  pressure  is  different  of  the  two  kinds  of  ports,  the  NINO  flow 


field  will  be  the  interdigitated  one.  From  the  analysis  above,  it 
is  found  that  the  parallel,  interdigitated  and  serpentine  flow  field 
can  all  be  gained  from  the  NINO  flow  field  at  some  extreme 
operation  condition.  Except  for  the  three  flow  fields,  there  are 
almost  infinity  flow  field  which  can  are  got  from  the  new  flow 
fields,  if  the  operation  condition  was  changed.  Thereby  which 
is  the  best  one  becomes  the  question  that  attracts  us  very  much. 
To  answer  this  question,  a  new  mathematic  model  is  presented 
as  follows. 

2.2.  Model  development 

The  PEMFC  model  described  here  is  developed  to  study  the 
effect  of  the  porous  layer  in  the  flow  field  and  the  velocity  varia¬ 
tion  at  the  port.  Consequently,  to  avoid  complexities,  we  decided 
to  develop  the  half-cell  model  for  the  cathode  side  of  the  fuel 
cell.  It  was  assumed  that  gas  in  the  flow  field  was  air  on  the 
cathode  side.  The  air  enters  the  domain  at  gas  inlet  in  the  bipolar 
plate.  The  flow  field  in  the  bipolar  plate  helps  in  the  distribution 
of  reactant  gas  onto  the  surface  of  the  electrode. 

2.2.1.  Problem  domain 

The  base  dimensions  for  channel  width,  shoulder  width  and 
the  whole  flow  field  are  shown  in  Fig.  3.  Since  it  is  very  complex 
to  optimize  the  novel  flow  field,  some  base  cases  are  discussed. 
The  main  content  will  be  concentrate  on  the  method  of  finding 
the  answer.  The  velocity  at  ports  1,  3,  5,  ports  2, 4,  ports  6,  8,  10, 
was  set  at  umi,  um 2,  wout2,  respectively,  and  the  pressure  at  ports 
7,  9  is  1.00  atm.  The  proportion  of  the  permeability  in  the  chan¬ 
nel  to  the  permeability  under  the  shoulder  will  varied  from  20: 1 
to  100: 1 .  The  permeability  of  gas  diffuser  under  shoulder  and  in 
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Fig.  3.  The  base  dimensions  for  channel  width,  shoulder  width  and  the  whole  flow  field. 


channel  varies  from  0.5  x  1.76  x  10-11  to  1.3  x  1.76  x  10-11, 
but  the  proportion  will  not  change.  The  effect  of  width  of  shoul¬ 
der  and  channel  is  studied  also. 

2.2.2.  Model  assumptions 

The  numerical  model  assumes  that 

1.  The  PEM  fuel  cell  operates  under  steady  state  condition  [4]. 

2.  There  are  no  temperature  changes  in  the  PEM  fuel  cell,  i.e. 
isothermal  operation  [4], 

3.  Laminar  flow  exists  everywhere  in  the  gas  channels  and  all 
the  transport  equations  were  formulated  for  laminar  behavior 

[7]. 

4.  Porous  media  exists  in  the  gas  diffusers  is  isotropic. 

5.  Gas  in  flow  field  is  ideal  gas  and  uncompress,  because  the 
pressure  drop  throughout  the  flow  field  is  very  small. 

6.  Water  produced  from  reaction  only  exist  in  gas,  in  this  paper 
we  do  not  consider  the  liquid  water  in  channel  [4], 

Actually,  assumptions  1-4  are  usually  considered  as  standard 
assumptions  in  many  papers  [10-15]. 


effect  of  gravity  for  fluid  flow  through  a  porous  medium  is, 


kp\  dP 
p  dx 

(7) 

kp\  dP 

P  3  V 

(8) 

where  kp  i  is  the  permeability  in  channel,  //  the  viscosity  of  the 
fluid  in  the  medium,  and  P  is  the  total  gas  pressure.  Here  we 
assume  that  the  channel  is  filled  with  porous  material  with  high 
permeability.  In  fact,  the  difference  of  permeability  in  channel 
and  under  shoulder  represents  the  distribution  ability  of  flow 
medium  in  them.  Darcy’s  law  is  related  to  the  Hagen-Poiseuille 
equation  for  laminar  flow  in  the  gas  channel.  The  kp\  an  //  are 
assumed  to  be  constant. 

When  Eqs.  (7)  and  (8)  are  substituted  into  Eq.  (5),  we  obtain, 
d(-(pkpl/p)(dP/  dx))  d(—(pkv\/p)(dP/dy)) 

-  =  Ochl 

dx  dy 

(9) 


2.2.3.  Governing  transport  equations 
The  reaction  equation  of  cathode  is 

O2  +  4H+  +  4e_  — »■  2H20  (4) 


The  conservation  of  mass,  also  called  equation  of  continuity 
in  channel  is  given  by  [4]: 


3 (pvx)  d(pvy)  _ 

—  *3ch  1 

dx  dy 


(5) 


where  p  is  the  density  of  the  fluid  in  the  medium,  vx  and  vy  are 
velocities  in  the  x-  and  y-direction  and  ,S'chi  is  the  source  term. 
This  source  term  includes  the  H+  and  water  transferring  from 
anode, 


^chl  — 


M  [  |+  i  MwXi 


(6) 


F  is  the  faraday  constant,  7.  the  electro-osmotic  drag  coefficient 
of  water,  and  i  is  the  current  density.  Darcy’s  law  neglecting  the 


There  are  three  species  that  are  oxygen,  water  and  nitrogen  in 
an  air  cathode.  Oxygen  is  fed  to  the  electrode  interface,  water  is 
generated  at  the  gas  diffusion  layer  (GDL)/electrode  interface. 
When  oxygen  reacts  at  the  GDL/electrode  interface,  the  Tafel 
type  equation  is  obtained  to  describe  the  current  distribution 
along  the  catalyst  layer. 


(10) 


where  I0  is  the  exchange  current  density  of  oxygen,  c0  the  mole 
concentration  of  oxygen  gas  in  catalyst  layer,  c0,r  the  reference 
oxygen  concentration,  and  k  is  the  transfer  coefficient. 

If  the  average  oxygen  concentration  in  channel  and  its  diffu¬ 
sion  layer  is  coa,  and  the  distribution  is  linear  shown  in  Fig.  4, 
the  relationship  of  coa,  c0,  c0i,  c0 2  is 


Q>a 


((c0  +  c0i)d2)/2  +  ((c0i  +  c02)dl)/2 


(11) 


dl  +d2 
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Fig.  4.  Oxygen  concentration  (a)  in  channel  and  its  diffusion  layer,  and  (b)  in 
diffusion  layer  under  shoulder. 


where  dl  is  the  height  of  channel,  and  d2  is  the  thickness  of  GDL. 
According  Fick’s  law,  oxygen  flux  in  channel  can  be  determined 


3cch 

dx 


Oxygen  flux  in  GDL  is 


/Vgdi,  =  Dgdl 


3cgdl 

dx 


(12) 


(13) 


If  there  is  no  chemical  reaction  in  the  flow  field,  the  mass  balance 
used  for  oxygen  and  water  vapor  can  be  described  as  follows, 


Dj 


.  3  2Cj 
dx 2 


+  o 


.  d2 


Cj 

dy2 


—  v 


dci 

dx 


-  Vv 


dcj 

dy 


=  0 


(14) 


where  i  represent  oxygen  or  water,  D'j  is  the  effective  molec¬ 
ular  diffusion  coefficient  of  species  i  in  the  gas  stream,  and  c; 
is  the  average  mol  concentration  of  species  i  in  channel.  The 
relationship  of  7)®  and  Dj  which  is  the  diffusion  coefficient  in  a 
nonporous  system  is  given: 

Of  =  Dte1-5  (15) 


where  e  is  the  porosity  of  the  channel.  In  this  model,  the  porosi¬ 
ties  and  the  effective  diffusion  coefficients  of  the  gases  are 
assumed  to  be  constant  and  invariant  of  pressure.  The  relation 
between  temperature  and  the  diffusion  coefficient  is  given, 

A  =  A',r(£)  (16) 

where  D,  r  is  the  reference  oxygen  or  water  diffusion  coefficient, 
T  the  operation  temperature,  and  7’,  is  the  reference  temperature. 


From  Fig.  3,  when  the  chemical  reaction  takes  place  at  any 
point  of  the  flow  field  plane,  the  oxygen  balance  is  got  from  Eq. 
(14): 


e  d"coa  e  3“Coa 
Ur*  UQ~ 


dx2 


dy2 


9c0a 

dx 


9c0; 
’  dy 


+  4  F  ° 


(17) 


where  U4F  represent  the  consumption  of  oxygen,  coa  represents 
mole  concentration  of  oxygen.  The  equations  of  water  vapor 
balance  is  got  from  Eq.  (14)  also: 


D[ 


9  iwa 

dx2 


D\ 


32Cwa 

'  dy2 


9cw 


3cw 


dx 


dy 


(1  +  27.); 
2  F 


=  0 
(18) 


where  (1  +  2X)i/2F  represents  the  production  of  water  vapor, 
and  cwa  represents  mole  concentration  of  water. 

Under  shoulder  the  conservation  of  mass  equation  is 


d(pvx)  d (pvy)  _  c 

„  „  —  ■Jshl 

dx  dy 


(19) 


5'ch  I 


(20) 


Darcy’s  law  neglecting  the  effect  of  gravity  for  fluid  flow 
through  a  porous  medium  is, 


vx 


Vy 


kp2  dP_ 
/x  dx 
k^d_P_ 
f  3 y 


(21) 

(22) 


The  kp2  and  //  are  assumed  to  be  constant.  When  Eqs.  (21)  and 
(22)  are  substituted  into  Eq.  (19),  we  obtain. 


d(-(pkp2/pi)(dP/dx))  d(—(pkp2/  p,)(dP/dy)) 

dx  dy 


Sshl 


(23) 


The  Tafel  type  equation  is  obtained  to  describe  the  current  dis¬ 
tribution  along  the  catalyst  layer. 


•  ,  c° 

;  =  I0 - exp 

co,r 


(24) 


If  the  average  oxygen  concentration  in  diffusion  layer  is  ca,  and 
the  distribution  is  linear,  the  relationship  of  coa,  c0  and  c0i  shown 
in  Fig.  4b  is 


Coa  — 


Co  +  Col 

2 


(25) 


According  Fick’s  law,  oxygen  flux  in  GDL  can  be  determined 

3cgdl 
dx 


Ngdl  —  Dgdi 


(26) 


The  equations  of  oxygen  and  water  balance  in  GDL  are  the  same 
as  the  ones  in  channel. 


2.2.4.  Boundary  conditions 

There  are  four  shoulders  that  can  be  changed  as  we  need  in 
this  flow  field.  Ports  1,  3,  5  having  the  same  flux  velocity  belong 
to  the  inlet  1;  ports  2,  4  belong  to  the  inlet2  which  have  two 
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times  flux  velocity  to  inletl;  ports  6,  8,  10  having  the  same  flux 
velocity  to  inlet2  belong  to  the  outlet  1;  ports  7,  9  belong  to  the 
outlet2  whose  pressure  was  set  to  1  atm.  All  ports  have  the  same 
width  LI  and  all  shoulders’  width  is  L2.  The  whole  length  and 
width  of  flow  field  is  L3  and  L4. 

In  order  to  make  the  result  figures  more  clear  and  reduce 
calculation,  symmetry  boundary  condition  was  applied  in  this 
paper,  and  part  of  the  flow  field  shown  in  Fig.  3  was  simulated. 


In  Fig.  3, 

AB  =  CD  =  Ll/2  (27) 

EF  =  GH  =  Ll/2  (28) 

BC  =  GF  =  L2  (29) 

DE  =  AH  =  L3  (30) 


At  the  inlet2  AB,  where  the  reactant  gases  are  fed  to  the  elec¬ 
trode,  the  pressure  and  the  mole  concentration  of  each  species 
are  specified  as  follows. 


Jn2 


in2 


dp 

dx 


kp\ 


uin2 


(31) 


At  the  inletl  CD,  the  mole  concentration  of  each  species  are 
specified  as  follows, 


in  I 


Jnl 


8P  (JL 

-  = - Ml 

dx  kpi 


(32) 


At  the  outlet  1  EF,  where  the  reactant  gases  flow  out  of  the  elec¬ 
trode,  the  pressure  and  the  mole  concentration  of  each  species 
are  specified  as  follows, 


dc0 

dx 


=  0, 


dcw 

dx 


=  0, 


P=Po 


(33) 


At  the  outlet2  GH,  the  mole  concentration  of  each  species  are 
specified  as  follows. 


dc0 

dx 


=  0, 


dcw 

dx 


=  0, 


dp 

dx 


F 

.  Mout2 
kpi 


(34) 


At  the  boundaries  AH,  DE,  Newman  boundaries  condition  was 
used 


dca 

dy 


=  0, 


dcw 

dv 


=  0, 


At  the  boundaries  BC,  FG 


dc0 

dx 


=  0, 


dev. 

dx 


=  0, 


(35) 


(36) 


where  M0,  Mw  and  Mn  represent  molecular  weight  of  oxygen, 
water  and  nitrogen,  respectively.  cna  is  the  average  nitrogen 
concentration.  Then  the  velocity  in  flow  field  will  be  got  from 
Darcy’s  law. 

Next  Eq.  (17)  of  oxygen  concentration  will  be  solved.  But 
the  current  density  i  is  a  function  of  oxygen  concentration,  and 
the  function  in  channel  is  different  from  the  one  under  shoulder. 

In  channel,  the  oxygen  flux  in  channel  and  in  GDL  are  equal 
to  each  other,  and  the  relationship  of  oxygen  concentration  in 
channel  and  in  GDL  can  be  got  from  Eqs.  (12)  and  (13), 


Ach  —  Dc  h 


9cch 

dx 


=  Ngdl  =  Dsh 


9cgdl 

dx 


(38) 


Assuming  the  concentration  variation  is  linear,  Eq.  (38)  can  be 
reduced 


_  Q>2  C01  _  C01  C0 

Dch - - -  =  DgdL- 


dl 

In  addition. 


d2 


N GDL  =  =  4ttt—  exP 

4  F  4Fcor 


'4kF 


v  RT 

Combining  Eqs.  (11),  (39)  and  (40) 
2 


Co  = 


.  Coa 


4 +  £1  +  3^2 
where  dimensionless  parameter 


c0dl  (  4k  F 

£1  =  4  77^ - exP  '  - 


4FDchc0,r  V  RT 


£2  =  4 


c0d2 


4/-/7gdi.c0, 


exp 


(4kF 


(39) 


(40) 


(41) 


(42) 


(43) 


Under  shoulder,  there  is  only  GDL.  Combining  Eqs.  (25),  (26) 
and  (40),  the  concentration  coa,  c0 


Co_  2  +  £2Coa 


(44) 


After  solving  Eq.  (17),  the  result  of  oxygen  concentration 
was  applied  to  solve  Eq.  (18)  to  get  the  distribution  of  water 
concentration. 

The  equations  were  solved  iteratively  with  finite-difference 
method,  and  the  solution  was  considered  to  be  convergent  when 
the  relative  error  between  two  consecutive  iterations  was  less 
than  10-5.  Stringent  numerical  tests  were  performed  to  ensure 
that  the  solutions  were  independent  of  the  grid  size.  An  80  x  60 
(length  x  width)  mesh  was  found  to  provide  sufficient  spatial 
resolution. 


3.  Numerical  techniques  and  procedures 

The  model  equations  developed  above  are  coupled  closely, 
so  the  whole  set  of  equations  must  be  solved  simultaneously  and 
iteratively.  First  we  solved  the  pressure  Eqs.  (9)  and  (23)  with 
the  concentration  of  water  and  oxygen  was  set.  The  gas  density 
can  be  calculated  by 

(37) 


4.  Results  and  discussion 

In  parallel  flow  field,  the  fluid  mainly  flows  in  channels,  and 
the  velocity  under  shoulders  is  nearly  zero.  Diffusion  is  the  main 
way  for  the  components  to  transport  between  the  place  in  chan¬ 
nel  and  under  shoulder.  So  reducing  the  width  of  shoulder  is 
advantageous  to  components’  transportation,  but  narrow  shoul¬ 
ders  will  increase  difficulty  of  manufacture. 


P  —  (44  —  Mw)c0  a  +  (44  —  Mw)cna  +  PMW/ RT 


1216 


X.  Zhou  et  al.  /  Journal  of  Power  Sources  158  (2006)  1209-1221 


In  interdigitated  flow  field,  the  fluid  flows  not  only  in  chan¬ 
nels  but  also  under  shoulders.  The  pressure  difference  between 
two  neighbor  channels  forces  the  fluid  to  flow  under  shoulder, 
and  the  components’  transportation  under  shoulder  in  interdig¬ 
itated  flow  field  becomes  much  better  than  that  in  parallel  flow 
field.  As  a  result,  the  interdigitated  flow  field  performs  better 
at  high  mass-transport  limiting  current  density  as  compared  to 
the  parallel  flow  field  [3],  But  the  velocity  at  the  end  part  of 
channels  (around  b  position  shown  in  Fig.  lb)  is  much  slower 
than  that  at  the  start  part  of  channels  (around  a  position  shown 
in  Fig.  lb),  and  the  reactant  concentration  at  this  part  will  be  rel¬ 
ative  low.  What’s  more,  high  velocity  in  GDL  and  catalyst  layer 
cause  the  catalyst  layer  and  membrane  dehydrate,  and  increase 
ohmic  overpotential,  when  no  water  is  added  to  humidify  the 
cathode  inlet  air.  So  the  interdigitated  flow  field  performs  worse 
at  low  or  middle  current  density  as  compared  to  the  parallel  flow 
field  [2],  In  addition,  in  order  to  keep  a  certain  fluid  flux,  higher 
pressure  must  be  applied  to  the  interdigitated  flow  field  with  the 
increasing  of  lost  power  consumed  by  fuel  cell  system  itself. 

Combining  the  two  flow  fields,  a  new  flow  field  is  shown  in 
Fig.  lc  forms.  Because  of  velocity  difference  of  the  neighbor 
inlets  and  outlets,  the  pressure  difference  between  the  neigh¬ 
bor  channels  is  generated  to  force  part  of  fluid  to  flow  through 
the  GDL  under  shoulder.  So  the  two  neighbor  channels  will 
served  as  interdigitated  flow  field,  and  the  shear  force  of  gas 
stream  helps  removing  the  liquid  water  result  in  reducing  the 
electrode-flooding  problem,  and  the  oxygen  concentration  will 
increase.  On  the  other  hand,  the  velocity  in  channels  is  rela¬ 
tive  high  because  of  higher  velocity  at  each  outlet  (end  part  of 
channels). 


Because  each  port  in  NINO  is  independent,  every  port  can  be 
fixed  to  be  inlet  or  outlet.  From  studying  of  NINO,  almost  infinity 
flow  field  modes  will  be  obtained,  and  through  the  studying  of 
these  modes,  it  is  forecast  that  some  outstanding  flow  field  will  be 
found  to  be  suit  to  PEMFC.  The  flow  field  with  two  kinds  of  inlets 
(inlet  1  and  inlet2)  and  two  kinds  of  outlets  (outlet  1  and  outlet2) 
was  studied  here.  The  parameters  were  shown  in  Table  1.  Inlet  1 
to  outlet  1  and  inlet2  to  outlet2  served  as  a  parallel  and  have  high 
velocity.  As  a  result,  the  oxygen  in  channel  renews  very  fast  and 
current  under  channel  will  increase.  Especially,  when  the  feed 
is  not  gas  but  methanol,  a  mount  of  CO2  gas  will  be  generated 
at  high  current  density,  and  the  channel  will  be  filled  with  CO2 
gas  [16].  At  this  time,  removing  CO2  gas  instead  of  increasing 
inlet  velocity  becomes  important.  Inlet  1  to  outlet  1  and  inlet2  to 
outlet2  can  achieve  this  goal,  because  relative  high  velocity  in 
channel  from  inlet  to  outlet. 

In  order  to  describe  the  flow  field  and  optimize  the  NINO 
flow  field,  a  new  simple  model  is  given.  Considering  the  com¬ 
plex  of  optimization,  some  simple  conditions  were  studied  and 
optimizing  work  will  be  reported  later. 

4.1.  Model  validation 

To  validate  the  current  two-dimension  model,  comparisons 
are  made  to  the  experimental  data  of  Ticianelli  et  al.  [17]  for  a 
single  cell.  The  schematic  of  the  computation  domains  is  pre¬ 
sented  in  Fig.  3,  and  the  basic  parameters  used  in  the  model 
validation  are  given  in  Table  1 . 

The  mathematic  model  developed  above  can  produce  cathode 
(or  anode  by  modification)  polarization  data  of  voltage  versus 


Table  1 

Conditions  and  parameters 

Parameter  Value  Source 


Number  of  the  shoulders 

Width  of  the  channel  beside  portl  and  port5,  LI  (m) 

Width  of  the  shoulder,  L2  (m) 

Length  of  the  whole  flow  field,  L3  (m) 

Width  of  the  whole  flow  field,  L4  (m) 

Width  of  the  portl,  L5  (m) 

Width  of  the  port5,  L6  (m) 

High  of  channel,  dl  (m) 

Thickness  of  GDL,  d2  (m) 

Air  viscosity,  //.  (N  m-2  s) 

Permeability  of  gas  diffuser  under  shoulder,  kp2  (m2) 

Permeability  of  gas  diffuser  in  channel,  kpl  (m2) 

Gas  diffuser  porosity,  e 

Oxygen  diffusivity  in  gas,  D0  (m2  s-1) 

Water  diffusivity  in  gas,  Z)w  (m2  s-1) 

Transfer  coefficient  of  oxygen  reduction  reaction,  k 
Reference  oxygen  concentration,  cQ,r  (molm-3) 

Exchange  current  density  of  oxygen  concentration  reaction,  IQ  (A  m-2) 
Faraday  constant,  F  (C  mol-1) 

Operating  temperature,  T  (K) 
b 

Cathode  overpotential,  r]  (V) 

Atmospheric  pressure,  Pq  (Pa) 

Cathode  oxidant  stoichiometric  flow  ratio 
Pressure  at  outlet  1  (Pa) 

Electro-osmotic  drag  coefficient  of  water 


4 


3  x  ltr3 

3  x  io-3 

27  x  nr3 

27  x  nr3 

i.5  x  ur3 

i.5  x  ur3 

1  x  It)-3 

3  x  10“4 

1.38  x  10-5f77298)102 

[22] 

1.76  x  10"11 

[4] 

1.76  x  1(T9 

Assuming 

0.4 

[23] 

1.775  x  10'5(r/273.15)15 

[4] 

2.56  x  10“5(r/307.15)15 

[4] 

0.5 

[4] 

34.5 

Calculated  value 

100 

[24] 

97487 

353.15 

1.5 

[7] 

0.1 

101325 

1.68 

Assuming 

Po 

1 

[3] 
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current  density.  To  validate  the  model  with  data  of  single  fuel 
cell,  the  other  parts  of  fuel  cell  were  assumed  to  be  a  resis¬ 
tance  because  of  the  high  activity  of  transport  and  reaction  for 
hydrogen  such  as  high  diffusion  coefficient  and  high  exchange 
current  [3].  The  internal  conductivity  of  the  fuel  cell  was  set  to 
be  4  x  104  S  m2  at  353  K  from  the  experiment  data  fitting,  and 
it  varied  with  temperature  by  [18]: 


a  =  exp 


j  (5139/1  —  3260) 


(45) 


where  /I  is  the  water  content,  and  it  is  set  to  be  constant.  Com¬ 
bining  the  conductivity  got  from  fitting  result  and  Eq.  (45),  the 
conductivity  at  different  temperature  is 


a  —  4  x  104  exp 


1268 


/  1 


\  353.1 5 


(46) 


Flow  field  length  /  mm 


xIO5 
1.014, 


S1 


Fig.  6.  Pressure  distribution  in  the  NINO  for  the  base  case. 


The  exchange  current  density  of  oxygen  reaction  is  given  by 
[19]: 

I0  =  l03.507-4001/r  (47) 

The  relationship  of  voltage  and  current  density  is 

E  =  E0-  rj c  -  ir  (48) 

where  E  is  the  voltage  of  cell,  Eq  the  thermodynamic  open  cir¬ 
cuit  potential,  i  the  average  current  density,  and  r  the  internal 
resistance  of  cell. 

The  curves  in  Fig.  5  are  the  numerically  prediction  showing 
good  agreement  with  the  newly  selected  experimental  voltage 
versus  current  density  at  the  different  operating  temperature, 
then  the  model  is  proved  to  be  an  accurate  model. 

4.2.  Base  case 

Fig.  6  shows  the  pressure  distribution  all  through  the  flow 
field.  The  pressure  difference  is  smooth  in  channel  and  sharp 
under  shoulder.  It  is  mainly  because  the  permeability  (kpl)  in 
channel  is  100  times  bigger  than  under  shoulders  (kp2)  and  pres¬ 
sure  difference  between  two  channels  is  obvious. 


Fig.  5.  Comparison  of  numerically  predicted  polarization  curves  with  experi¬ 
mental  data. 


The  profiles  of  x-directions,  y-d  i  recti  on  s  and  average  veloc¬ 
ity  are  given  in  Figs.  7-9.  The  x-direction  velocity  in  channel 
decreases  along  x-direction  in  high-pressure  channel  and  it  is 
much  higher  than  that  under  shoulders.  The  highest  y-direction 
velocity  appears  in  channel  near  the  ports,  because  the  flow  rate 
near  the  ports  is  much  bigger  than  that  at  other  parts  of  the  flow 
field.  Here  the  pressure  difference  in  y-direction  under  shoul¬ 
ders  is  sharp.  As  a  result,  the  velocity  under  shoulder  is  mainly 
along  y-direction.  When  the  velocity  of  the  ports  varies,  the  inlets 


Fig.  7.  x-Direction  velocity  distribution  in  the  NINO  for  the  base  case. 


Fig.  8.  y-Direction  velocity  distribution  in  the  NINO  for  the  base  case. 
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and  outlets  will  transform  and  the  velocity  under  shoulder  will 
change.  Using  this  method  all  kinds  of  different  flow  fields  can 
be  studied,  and  some  good  ones  will  be  selected. 

The  oxygen  distribution  is  shown  in  Fig.  10.  The  concentra¬ 
tion  of  oxygen  decreases  along  the  channel  and  reach  the  lowest 


Fig.  10.  Oxygen  distribution  in  the  NINO  for  the  base  case. 


near  the  outlet2.  Since  the  velocity  in  channel  is  very  high,  the 
concentration  decreases  slowly.  But  the  concentration  changes 
very  fast  under  shoulders,  because  of  the  consumption  and  slow 
velocity.  In  addition,  we  know  from  Eq.  (10)  that  the  current 
distribution  corresponds  to  the  oxygen  distribution  shown  in 
Fig.  10.  The  shape  of  this  figure  is  similar  to  the  result  of  con¬ 
ventional  flow  field  [20].  But  the  vale  of  oxygen  concentration 
under  shoulder  is  not  symmetric  because  of  y-direction  velocity. 

The  water  vapor  distribution  shown  in  Fig.  1 1  has  the  oppo¬ 
site  trends.  The  water  vapor  generated  in  flow  field  comes  from 
two  ways:  first,  the  water  was  generated  at  the  place  where  the 
oxygen  was  consumed;  second,  net  water  transports  with  proton 
transportation  from  anode  to  cathode.  The  velocity  under  shoul¬ 
ders  is  slower  than  in  channels,  and  the  water  concentration 
under  shoulders  becomes  higher  than  in  channel.  The  highest 
concentration  of  water  appears  near  the  outlet2.  When  the  elec¬ 
trode  overpotential  is  changed,  and  the  corresponding  current 
densities  are  obtained  from  the  following  equation: 


i(x,  y)  ds 


(49) 


where  i  represents  the  average  current  of  the  whole  fuel  cell, 
A  represents  the  area  of  the  electrode,  and  i(x,  y)  is  the  current 
function.  The  polarization  curve  for  the  base  case  shown  in  Fig.  5 
is  also  obtained.  The  average  current  density  for  the  base  case 
is  0.9226  A  cm-2. 


4.3.  Effect  of  permeability  of  gas  diffuser  (kp2)  under 
shoulder 

Fig.  12  shows  the  effect  of  kp2  variation  on  the  pres¬ 
sure  distribution  at  x  =  0.5L3.  The  values  of  kp2  include 
0.5  x  1.76  x  10“n,  0.7  x  1.76  x  10“n,  0.9  x  1.76  x  10“u, 
1.1  x  1.76  x  10-11  and  1.3  x  1.76  x  10-11  with  the  other 
parameters  keeping  as  base  case.  The  proportion  of  kpl  :kp2  will 
not  be  changed.  The  distributions  of  velocity,  oxygen  concen¬ 
tration  and  water  concentration  have  the  same  value  calculated 
from  the  model,  but  the  pressure  decreases  dramatically  in  high- 
pressure  channel  shown  in  Fig.  12  when  kp2  increased.  This  is 


4  3 

Width  /  mm 


Fig.  11.  Water  vapor  distribution  in  the  NINO  for  the  base  case. 


X.  Zhou  et  al.  /  Journal  of  Power  Sources  158  (2006)  1209-1221 


1219 


0.315 


- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

0  1  2  3  4  5  6 

Flow  Field  width  /  mm 

Fig.  12.  Effect  of  permeability  of  gas  diffuser  under  shoulder  on  pressure  dis¬ 
tribution  at  x  =  0.5L3. 

mainly  because  of  applying  of  the  velocity  boundary  conditions. 
The  dramatic  variation  of  pressure  can  be  easily  explained  by 
Eqs.  (7)  and  (8). 

4.4.  Effect  of  the  proportion  ofkpl  in  channels  to  kp2 
under  shoulders  (kpl:kp2) 

The  effect  of  kpl:kp2  (kpl  represent  the  kp  in  channels, 
kp2  represent  the  kp  under  the  shoulders)  is  investigated.  When 
kpl:kp2  is  in  case  of  20:1,  40:1,  60:1,  80:1,  100:1,  the  effect 
of  these  variations  was  studied  with  the  other  parameters  being 
same  as  base  case.  The  effects  of  vary  aspects  are  shown  from 
Figs.  13-16.  When  we  changed  the  proportion,  the  kp2  was  kept 
in  1.76  x  10“ 1 1  m2.  Fig.  13  shows  that,  as  the  increasing  of  the 
proportion,  the  pressure  increases  both  in  high-pressure  chan¬ 
nel  and  in  low-pressure  channel.  Especially  at  the  proportion  of 
20:1,  the  pressure  difference  trends  to  be  small.  That  is  to  say, 
increasing  the  proportion  can  undulate  the  pressure  distribution 
in  the  NINO  flow  field.  But  this  trend  becomes  weak  when  the 
proportion  of  kpl  to  kp2  increases  up  to  60: 1 .  The  effect  on  the 
velocity  was  shown  in  Fig.  14,  the  velocity  increases  with  the 
rising  ofkpl  in  channel  while  decreases  under  shoulder,  and 


Fig.  13.  Effect  of  the  proportion  of  kpl  to  kp2  on  pressure  distribution  at 
je  =  0.5L3. 
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Fig.  14.  Effect  of  the  proportion  of  kpl  to  kp2  on  average  velocity  distribution 
at  x  =  0.5L3. 


Width  /  nm 

Fig.  15.  Effect  of  the  proportion  of  kpl  to  kp2  on  concentration  of  oxygen 
distribution  at  x  =  0.5L3. 

the  quantity  of  the  increasing  becomes  large  near  20:1.  It  can 
be  seen  from  Fig.  13  that  the  pressure  difference  between  two 
neighbor  channels  increases  as  increasing  kpl.  But  the  velocity 
under  the  shoulders  decreases  with  the  rising  of  kp  1  because  the 
velocity  in  x-direction  increases  faster. 


Width  /  mm 

Fig.  16.  Effect  of  the  proportion  ofkpl  to  kp2  on  concentration  of  water  distri¬ 
bution  at  v- 0.5L3. 
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The  concentrations  of  oxygen  and  the  water  were  shown  in 
Figs.  15  and  16.  However  the  distributions  of  oxygen  and  water 
also  have  relationship  with  the  velocity  at  the  ports.  The  velocity 
we  selected  here  is  different  in  neighbor  channels  and  can  help 
to  distribute  the  oxygen  in  the  NINO  flow  field.  When  the  kpl 
increases,  the  velocities  decrease  in  channel  and  increase  under 
shoulder.  However  the  state  of  water  is  reverse,  because  of  the 
constantly  producing  water. 

The  polarization  curve  under  the  five  conditions  was  not  given 
here,  because  the  oxygen  concentration  in  each  condition  is 
almost  same  to  each  other.  It  is  mainly  because  the  decreasing  of 
oxygen  concentration  under  shoulder  counteracts  the  increasing 
of  it  in  channel  with  the  decreasing  of  kpl:kp2. 

Therefore  it  is  helpful  to  enhance  the  PEMFC  performance 
in  moving  the  water,  as  increasing  the  proportion  of  kpl  to 
kp2. 

4.5.  Effect  of  width  of  channel  and  shoulder 

Effect  of  width  of  shoulder  and  channel  was  shown  in 
Figs.  17-21.  The  width  summation  of  the  neighbor  channel 
and  shoulder  was  kept  in  constant  value;  the  proportion  of  the 
channel  and  shoulder  ( Wcha'.  VFsh0)  was  changed  with  the  other 
parameters  keeping  as  base  case.  The  pressure  distribution  of 
the  NINO  flow  field  at  a  =  0.5L3  is  shown  in  Fig.  17.  When 
the  width  of  channel  decreases,  the  pressure  increases  in  high- 
pressure  channel,  while  decreases  in  low-pressure  channel  until 
Wcha:Wsho  is  1:1  and  then  increases. 

In  addition,  we  can  mention  that  the  pressure  difference  varies 
with  the  changing  of  proportion.  This  causes  the  velocity  under 
all  shoulders  shown  in  Fig.  18  increase  until  1:2.  However  the 
velocity  in  channels  increases  all  the  time.  Because  of  the  varia¬ 
tion  of  the  velocity,  the  distributions  of  oxygen  and  water  shown 
in  Figs.  19  and  20  must  redistribute.  The  concentration  of  oxy¬ 
gen  increases  only  in  high-pressure  channel,  but  decreases  in 
low-pressure  channel  and  under  shoulder.  When  the  width  of 
channel  increases,  the  fluid  in  the  NINO  flow  field  becomes 
homogeneous.  However,  the  homogeneous  distribution  of  water 
has  the  reverse  trends. 


Fig.  17.  Effect  of  width  of  shoulder  and  channel  on  the  pressure  distribution  at 
x  =  0.5L3. 


Fig.  18.  Effect  of  width  of  shoulder  and  channel  on  the  velocity  distribution  at 
*  =  0.5L3. 

The  polarization  curve  was  shown  in  Fig.  21.  The  con¬ 
tinually  increasing  of  the  current  has  direct  relation  with  the 
decreasing  of  width  of  channel.  The  intrinsic  reason  should 
be  the  high  mass-transport  activity  under  shoulder,  though  the 
concentration  under  shoulder  decreases  with  the  decreasing 
of  Wcha  ■  Wsho  -  This  result  is  not  like  the  result  in  the  three- 
dimension  model  where  there  is  an  optimal  proportion  of  width 
of  channel  to  the  width  of  shoulder  [21].  It  is  mainly  because 
of  the  assumption  that  the  channel  and  GDL  have  the  same 
height. 


Fig.  19.  Effect  of  width  of  shoulder  and  channel  on  the  concentration  of  oxygen 
atx  =  0.5L3. 
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Fig.  20.  Effect  of  width  of  shoulder  and  channel  on  the  concentration  of  water 
at.r  =  0.5L3. 


Current  density  /  Anv2 

Fig.  21.  Effect  of  width  of  shoulder  and  channel  on  polarization  curve  of  a 
PEMFC. 

5.  Conclusions 

The  new  flow  field  presented  here  can  easily  be  transformed 
to  the  general  one  such  as  parallel  and  interdigitated  flow  field. 
There  are  2 n  independent  ports  in  the  NINO.  The  ports  in  this 
paper  are  divided  into  four  kinds:  inlet  1,  inlet2,  outletl,  outlet2 
by  different  boundary  condition  at  the  ports.  A  two-dimension 
mathematic  model  was  developed  to  describe  the  distribution 
of  pressure,  x-di  recti  on  velocity,  y-di  recti  on  velocity,  average 
velocity,  oxygen,  water  and  the  polarization  curve  of  the  air 
cathode,  which  constructed  by  the  NINO. 

The  numerically  prediction  shows  good  agreement  with  the 
newly  selected  experimental  voltage  versus  current  density  at 
the  different  operating  temperature. 

With  the  increasing  of  kp,  the  pressure  decreases  dramatically 
in  high-pressure  channel.  But  the  other  results  such  as  x-direction 


velocity,  y-direction  velocity,  average  velocity,  oxygen,  water 
and  the  polarization  curve  will  not  change. 

When  the  proportion  of  kpl  to  kp2  increases,  the  velocities 
decrease  in  channel  and  increase  under  shoulder.  As  a  result, 
the  refreshment  of  oxygen  increases  fast  in  high-pressure  chan¬ 
nel  and  decreases  under  shoulder  and  in  right  channel,  and  the 
state  of  water  is  reverse.  But  it  scarcely  affects  the  polarization 
curve. 

With  the  decreasing  of  VFcha:VFSh0,  the  current  continually 
increases  because  of  the  high  mass-transport  activity  under 
shoulder,  though  the  concentration  under  shoulder  decreases. 
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